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An updated next-to-leading order (NLO) QCD analysis of all presently available longitudinally polarized deep- 
inelastic scattering (DIS) data is presented in the framework of the radiative parton model. 



1. INTRODUCTION 

Many new or improved measurements of the 
spin asymmetry 

9i(x,Q 2 ) 



F 2 (x,Q 2 )/[2x(l + R(x,Q 2 ))Y 



(1) 



R = Fl/2xFi, in longitudinally polarized DIS 
off various different targets (p, n, d) have become 
available in the past two years: El 42, El 43, and 
SMC have presented their final data sets [ [j] , and 
new results from HERMES, E154, and E155 have 
been reported recently [ ||. Since none of these 
measurements were included in our original NLO 
QCD analysis [ |j, it seems to be appropriate to 
reanalyse these data^j within the framework of [ 

In NLO the polarized structure function g\ in 
(0) reads (suppressing all x and Q 2 dependence) 
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(Aq + Aq)® (l + ^AC, 
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where AC 9iff are the spin-dependent Wilson coef- 
ficients and the symbol ® denotes the usual con- 
volution in x space. From (^]) it is obvious that 
inclusive DIS data [ Q, § can reveal only informa- 
tion on Aq + Aq but neither on Aq and Aq nor 
on A<?, which enters (Q) only as an 0(a s ) correc- 
tion. Thus one can either stick to a comprehen- 
sive analysis of polarized DIS [ f|, |) or one has to 
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impose certain assumptions about the flavor de- 
composition in order to be able to estimate pro- 
cesses other than DIS for upcoming experiments 
like RHIC. As in [ [| , and other more recent QCD 
analyses [ [| fjj, we follow the latter option. 

2. DETAILS OF THE GRSV ANALYSIS 

Due to the lack of space we shall be rather 
brief and concentrate only on the most impor- 
tant changes since [ The limited amount of 
data demands a reasonably simple, but flexible 
enough ansatz for the polarized densities such as 



Af(x,n 2 ) = N f x a f(l-xfff(x,» 2 ) 



(3) 



with / = u,u,d,d,s,s,g. Note that in [ ^| we 
actually used / = u v , d v in (^) instead of u and 
d, however the positivity constraint 



\Af(x,Q 2 )\<f(x,Q 2 ) 



(4) 



which we want to exploit in our analysis, does 
not necessarily hold for the valence densities. Of 
course, the bound ([|) is strictly valid only in LO 
and is subject to NLO corrections [ || because 
the Af become unphysical, scheme-dependent 
objects in NLO. However the corrections are not 
very pronounced, in particular at large x [ 0], the 
only region where (Q) imposes some restrictions 
in practice. We therefore use (Q) also in NLO. 

To further simplify (|^) we assume that Aq = 
Au = Ad, fix N u ,d by the relations between 
the first moments of the non-singlet combinations 
Aq3.8 and the F and D values (using the updated 
value for \g A /g v \ = 1.2670 ± 0.0035 [ §), and 
take As = As = XAq. In the latter relation we 
choose A = 1 (SU(3)f symmetric sea), but simi- 
larly agreeable fits are obtained, e.g., for A = 1/2 
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(see also [ 0), as well as by using an independent 
x shape for As, reflecting the above mentioned 
uncertainty in the flavor separation. 

For the unpolarized distributions / in and 
(^) we use the updated GRV densities [ [l(J and 
also adopt their values for the input scale /i (~ 
0.6 GeV in NLO) and a s {M 2 z ) = 0.114. Note that 
in [ [l^] the RG equation for a s is now solved 
exactly instead of using the approximative NLO 
formula (see, e.g., [[|), which is more appropriate 
at low Q 2 where many of the polarized data lie. 

To fix the remaining parameters in (^|) we per- 
form fits to the directly measured spin asymme- 
try (Q) in LO and NLO, which is mandatory if 
one wants to adopt these distributions in a con- 
sistent analysis of the perturbative stability of po- 
larized processes. Also most MC programs only 
contain LO matrix elements, and hence LO densi- 
ties are more appropriate. Apart from the above 
outlined 'standard scenario' fit, equally good fits 
can be performed in a 'valence scenario' (see 
which is based on the assumption [ [ll[] that the 
F, D values fix only the valence parts of Aq^^g. 
The main feature of such a model is the possibil- 
ity to describe the data with a vanishing strange 
sea. Due to the lack of space we do not pursue 
this scenario here and restrict ourselves in what 
follows to the results obtained in the NLO (MS) 
'standard scenario' framework. 

3. RESULTS 

A comparison of our new NLO fit with the 
available A\(x, Q 2 ) data [ 111 || and with our pre- 
vious analysis [E| is presented in Fig. 1. The total 
X 2 values of the new and old fits are 147.4 and 
183, respectively, for 185 data points and adding 
statistical and systematical errors in quadrature. 
As can be seen, sizeable differences appear only 
in case of the neutron asymmetry, which leaves 
its footprint also in the individual parton densi- 
ties A/ shown in Fig. 2. Since the neutron data 
mainly probe Ad, the most prominent changes 
are observed here. 

The differences in the sea and in Ag only reflect 
the fact that they are constrained to a much lesser 
extent by the data than u and d. In particular Ag 
remains to be hardly constrained at all, which is 
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Figure 1. Comparison of our reanalysed NLO 
results for A±(x, Q 2 ) (solid lines) with available 
data sets (not shown, but included in the fit: 
E155 (deuteron) and the old EMC and SLAC 
data) and with our previous analysis [ |i| (dashed 
lines) . 

not surprising due to the lack of any direct infor- 
mation on Ag so far. Therefore we show in Fig. 2 
also the results of two other fits, which are based 
on additional constraints on Ag. For the 'Ag = 0' 
fit we start from a vanishing gluon input in (0), 
and the 'static Ag' is chosen in such a way that 
its first moment becomes independent of Q 2 . This 
can be achieved by setting dAg(Q 2 )/d\nQ 2 = 
and yields in LO Ag ~ -AS/2 (see [ |||), where 
AE is the total helicity carried by quarks (the re- 
lation is only subject to a small NLO correction). 
Both gluons give also excellent fits to the avail- 
able data and do not affect the results for u and d 
(see Fig. 2). In fact we can obtain fits without any 
significant change in \ 2 f° r Ag([i 2 ) in the range 
-0.3. . . 0.6 (our best fit has 0.28, corresponding to 
Ag(10GeV 2 ) ~ 0.9). The uncertainty in Ag is 
compatible to the one found in [ |jj , although our 
gluons extend also to slightly negative first mo- 
ments. It is interesting to observe that for our 
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Figure 2. The polarized NLO MS densities at 
Q 2 — 4 GeV 2 as obtained in the new and old [ | 
GRSV analyses. Also shown are the distributions 
obtained in two other fits employing additional 
constraints on Ag (see text). In all fits shown we 
have chosen a flavor symmetric sea. 



best fit gluon the spin of the nuclcon 

Ag(Q 2 ) + L Z (Q 2 ) 



(5) 



is dominantly carried by quarks and gluons at our 
input scale /i, and only during the Q 2 evolution a 
large negative L Z (Q 2 ) is being built up in order 
to compensate for the strong rise of Ag(Q 2 ), see 
Fig. 5 in [ [T^]. For the 'static Ag' the situation 
is completely different: by construction the quark 
and gluon contributions to (^) cancel each other 
implying that for all values of Q 2 the spin is en- 
tirely of angular momentum origin, contrary to 
what is intuitively expected. 

Inevitably the large uncertainty in Ag implies 
that the small x behaviour of g\ is completely 
uncertain and not predictable (see Fig. 3 in [ [l2| ) , 
which translates also into a large theoretical error 
from the x — > extrapolation when calculating 
first moments of <?i . Taking our best fit we obtain 
= 0.133 and r™ = -0.624 at Q 2 = 10 GeV 2 in 
agreement with a recent SMC QCD analysis [ ||. 

The two challenging questions concerning po- 
larized parton densities are still Ag and the flavor 
decomposition. Recent semi-inclusive DIS results 
[ |f3) may help to unravel the latter, but major 
progress, in particular on Ag, can be only ex- 
pected from RHIC. A realization of the optional 



upgrade of HERA to a polarized collider would be 
also very helpful to gain more insight into the spin 
structure of the nucleon and the photon, with the 
latter being completely unmeasured so far. 
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